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1  Introduction 


We  describe  two  types  of  molecular  nitrogen  fluorescence  lidars  which  can  be  used  to  probe 
the  D  region  of  the  ionosphere  during  periods  of  moderate  or  strong  aurora]  activity.  Both 
types  of  lidars  measure,  as  a  function  of  altitude,  the  density  of  a  particular  excited  or 
ionized  state  of  Nj.  Excited  and  ionized  states  of  Nj  exist  in  aurora  due  to  impact  of 
ground  state  Nj  with  high  energy  electrons. 

The  lidars  can  also  be  used  for  high  spatial  resolution  diagnostics  during  ionospheric 
high  power  radio  frequency  heating  experiments.  In  this  application  both  lidars  utilize  the 
fact  that  ionospheric  heating  experiments  are  often  performed  during  periods  of  moderate 
auroral  activity.  Moderate  auroral  activity  is  necessary  for  there  to  be  generation  of  ELF 
and  VLF  radiation  by  ionospheric  electrons  during  rf  heating  [l].  During  the  application 
of  high  power  rf  the  temperature  of  the  electrons  increases  which  in  turn  causes  changes 
in  the  densities  of  excited  and  ionized  Nj  states.  The  lidars  then  detect  the  changes  in 
density  of  a  particular  excited  or  ionized  state.  Changes  in  electron  temperature  are  then 
inferred  from  these  changes  in  density. 

The  first  type  of  lidar  we  describe  is  a  conventional  laser  induced  fluorescence  (LIF)  li- 
dar  operating  in  the  1*‘  positive  system  (Nj(AsE^)-+Nj(BsIIa))  or  the  1**  negative  system 
(Nj  (X2Ea  )-»Nj  (J3lE+))  of  molecular  nitrogen  (see  the  enery  level  diagram  for  Nj  in  Fig¬ 
ure  1).  Several  different  excitation/fluorescence  wavelength  combinations  are  considered, 
depending  on  the  choice  of  laser. 

The  second  type  of  lidar  we  refer  to  as  a  pump/LIF  lidar.  Two  lasers  are  used.  The  first 
laser,  the  so-called  pump  laser,  pumps  the  excited  nitrogen  from  N2(A3E+)  to  N2(B3na) 
(1**  positive  system).  The  second  laser,  the  so-called  probe  laser,  is  then  used  to  do  LIF 
lidar  on  the  2nd  positive  system  (N2(Bsna)-+Na(C3no)).  As  a  comparison  we  show  that  a 
conventional  LIF  lidar  operating  in  the  2nd  positive  system  without  a  1“  positive  pump  is 
not  viable  for  D  region  studies. 

We  consider  several  different  lasers  for  both  types  of  lidars.  Figure  2  is  a  schematic 
of  an  energy  level  diagram  showing  four  different  lidar  scenarios.  We  will  consider  these 
four  scenarios  throughout  this  report.  The  first  three  scenarios  refer  to  LIF  lidars  and  the 
fourth  refers  to  a  pump/LIF  lidar.  For  each  scenario  is  list  of  lasers  we  consider  and  their 
nominal  wavelengths. 

All  lasers  of  both  lidars  can  be  nitrogen  lasers  operating  in  the  band  of  interest  [2]. 
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INTERNUCLEAR  DISTANCE  (A) 


Figure  1:  Energy  levels  of  molecular  nitrogen  [12].  The  electronic  states  of  interest  in  this 
report  are  highlighted. 
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Figure  2:  Four  lidar  scenarios  (first  three:  LIF;  fourth:  pump/LIF),  the  lasers  considered 
for  each,  and  the  nominal  wavelengths  at  which  they  operate. 
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This  precludes  the  need  for  wavelength  tuning.  Alternatively,  for  the  1*‘  positive  system 
we  consider  tunable  dye  and  Ti:Sapphire  lasers,  and  nontunable  Nd:YAG  and  Nd:YLF 
lasers.  The  latter  two  lasers  happen  to  operate  at  wavelengths  which  coincide  with  specific 
rotational  transitions  of  the  1'*  positive  system  (see  Figure  2).  Since  the  Phillips  Labo¬ 
ratory/Ionospheric  Effects  division  (PL/GPIM)  mobile  lidar  is  equipped  with  a  NdrYAG, 
this  type  of  laser  is  an  especially  attractive  choice.  For  the  1*  negative  and  2nd  positive 
systems  we  consider  dye  and  Ti:Sapphire  lasers. 

In  addition  to  detecting  fluorescence  both  types  of  lidars  can  simultaneously  detect 
elastic  (Rayleigh)  backscattered  laser  radiation  from  air  molecules  as  a  function  of  altitude 
up  to  and  including  the  D  region  [3].  This  provides  density  measurements  of  the  ambient 
atmosphere.  Then,  with  standard  lidar  techniques,  the  temperature  of  the  ambient  atmo¬ 
sphere  can  be  deduced  from  these  density  measurements.  The  elastic  lidar  can  be  used  to 
monitor  gravity  wave  activity  at  D  region  altitudes  which  may  have  significant  effects  on 
ionospheric  heating  [4]. 

Due  to  the  relatively  small  beam  divergence  of  a  lidar,  high  horizontal  spatial  resolution 
can  be  achieved  as  well  with  a  scannabie  lidar  system. 

We  discuss  both  types  of  lidars  in  detail  in  the  following  sections.  We  present  signal 
to  noise  calculations  and  we  discuss  the  reasons  for  choosing  the  various  transitions  and 
the  relative  merits  of  each  type  of  lidar.  The  results  of  our  signal  to  noise  calculations  are 
compiled  in  Table  8  on  page  19. 

Sections  2  and  3  discuss  LIF  lidar  and  pump/LIF  lidar  respectively.  Both  sections 
present  a  general  description  of  each  lidar  and  show  how  number  of  detected  photons  is 
calculated. 

Section  4  shows  how  the  density  of  the  aurorally  excited  or  ionized  molecular  nitrogen 
is  calculated.  These  densities  are  required  for  determining  number  of  signal  photons.  The 
density  calculations  are  based  on  observed  column  emission  data  from  aurora. 

Section  5  discusses  signal  to  noise  ratio  and  presents  results  of  signal  to  noise  calcu¬ 
lations.  Section  5.1  presents  the  number  of  photons  expected  for  each  lidar  scenario  that 
we  consider.  All  the  necessary  lidar  parameters  required  for  this  calculation  are  listed. 
Section  5.2  discusses  and  quantifies  the  noise  processes  that  we  consider.  Section  5.3  then 
presents  signal  to  noise  ratios  for  each  lidar  scenario. 

Appendix  Appendix  A  shows  how  we  calculate  wavelengths  of  individual  rotational 
lines  of  N2  and  N£. 
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Appendix  Appendix  B  summarizes  laboratory  based  laser  induced  fluorescence  exper¬ 
iments  performed  at  the  UCLA  Plasma  Physics  Laboratory  in  February,  1993. 

Appendix  Appendix  C  shows  passive  emission  data  acquired  with  the  receiver  of  PL/GPIM 
mobile  lidar  system  at  the  HIP  AS  observatory  in  Alaska. 

2  LIF  Lidar 

2.1  General  Description 

The  LIF  lidars  that  we  consider  use  a  narrow  band  laser  to  excite  D  region  molecular 
nitrogen  from  a  single  rotational  level  of  a  specific  vibrational  level  of  a  specific  electronic 
state  to  another  rotational  level  of  a  specific  vibrational  level  of  a  higher  energy  electronic 
state.  The  molecule  then  decays  (fluoresces)  to  a  rotational  and  vibrational  level  of  the 
original  electronic  state.  The  lower  rotational  and  vibrational  levels  may  or  may  not  be 
identical  to  the  original  ones.  Also,  the  lower  electronic  state  is  an  excited  state  of  N2  or  a 
state  of  N2.  A  telescope  situated  adjacent  to  the  laser  detects  the  fluorescence  radiation. 

In  designing  a  LIF  lidar  consideration  must  be  given  to  the  electronic  system  (i.e, 
electronic  transition),  the  excitation  and  fluorescence  bands  within  the  electronic  system 
(i.e.,  vibrational  transitions),  and  the  excitation  and  fluorescence  rotational  transitions. 

Electronic  system 

The  choice  of  electronic  system  is  driven  by  the  need  for  a  practical  wavelength,  a 
strongly  absorbing  and  radiating  system,  and  a  relatively  high  density  of  the  lower  elec¬ 
tronic  state.  We  consider  LIF  lidars  operating  in  the  1**  positive,  1*‘  negative,  and  2nd 
positive  systems.  The  1**  positive  system,  which  has  wavelengths  from  the  visible  red  to 
the  mid  infrared,  is  considered  because  of  the  relatively  high  number  density  of  N2(A3E+) 
in  the  D  region  during  an  aurora.  The  density  of  N2(A3E+)  is  high  because  of  its  relatively 
long  lifetime  at  D  region  altitudes.  For  example,  at  100  km  the  lifetime  of  N2(A3E+)  is 
approximately  45  msec.  The  dominant  loss  mechanisms  are  due  to  collisional  quenching 
with  O  and  02.  (As  a  comparison,  the  radiative  lifetime  of  N2(A3E+)  is  2  sec.)  The  density 
of  N2(A3E+)  is  derived  in  Section  4. 

We  consider  the  1**  negative  system,  which  has  wavelengths  from  the  near  UV  to  the 
blue  visible,  because  of  the  relatively  high  absorption  cross  sections  and  Franck-Condon 
factors  of  some  bands  in  this  system.  Also,  the  density  of  N2  (XxE+)  is  relatively  high, 
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although  not  as  high  as  the  density  of  N2(A3E+). 

A  LIF  lidar  operating  in  the  2nd  positive  system,  which  has  wavelengths  from  the  far 
UV  to  near  UV,  is  not  viable  for  D  region  studies  because  of  the  very  low  number  density 
of  N,(£SII,).  The  density  of  Nj(B3II,)  is  low  because  of  its  relatively  short  radiative 
lifetime  of  approximately  10  n sec.  However,  because  of  the  high  density  of  Na(AsE+)  the 
2nd  positive  system  can  be  used  in  LIF  lidar  if  some  N2(A3E+)  molecules  are  first  excited 
to  Na^II,).  This  is  the  basic  idea  of  a  pump/LIF  lidar  discussed  in  the  next  section. 
Here  we  present  conventional  2nd  positive  LIF  lidar  as  a  comparison. 


Lower  vibrational  and  rotational  levels  of  excitation  transition 


The  lower  vibrational  and  rotational  states  of  the  excitation  transition  are  chosen, 
where  possible,  to  be  those  which  are  most  abundant  at  D  region  altitudes  during  aurora. 
For  the  situations  in  which  the  lasers  are  not  tunable  (i.e.,  Nd:YAG  and  NdrYLF)  the 
vibrational  and  rotational  levels  are  pre-determined  by  the  laser  wavelength.  For  the 
situations  in  which  a  nitrogen  laser  is  used  there  may  be  several  lower  rotational  levels 
excited,  depending  on  the  laser  output.  In  some  situations  the  most  abundant  vibrational 
state  may  not  be  well  coupled  to  an  upper  state  (i.e.,  low  absorption  cross  section)  and  an 
alternative  lower  vibrational  state  must  be  considered. 

The  density  of  a  vibrational  state  of  an  electronic  system  is  determined  from  past 
observations  of  aurora  reported  in  the  literature  [5].  Section  4  describes  how  we  calculate 
the  density  from  these  observations. 

The  density  of  a  rotational  level  is  determined  by  assuming  the  rotational  states  of  a 
given  vibrational  level  to  be  in  thermal  equilibrium  with  the  surroundings.  In  this  case 
the  fraction  of  molecules  J2ei„(J)  in  electronic  state  e  and  vibrational  state  v  that  are  in 
rotational  state  J  at  temperature  T  is  [6] 


R...V) 


hcBv 

kT 


(2  J  4- 1)  exp 


hcBvJ[J  ■+■ 1) 
kT 


(1) 


where  the  constant  Bv  is  a  function  of  spectroscopic  constants  and  the  vibrational  quantum 
number  (see  Appendix  Appendix  A). 


Upper  vibrational  and  rotational  levels  and 

lower  vibrational  and  rotational  levels  of  the  fluorescence  transition 

Ideally,  the  upper  vibrational  level  and  the  lower  vibrational  level  of  the  fluorescence 
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transtion  are  chosen  to  be  those  levels  which  give  the  greatest  product  of  absorption 
cross  section  and  Franck-Condon  factor.  The  former  is  proportional  to  the  strength  of 
the  excitation  transition  and  the  latter  is  proportional  to  the  strength  of  the  fluorescence 
transition.  In  reality,  however,  consideration  must  be  given  to  the  wavelengths  of  both 
transitions  and  to  the  lifetime  of  the  upper  state.  The  choice  of  excitation  wavelength  is 
driven  by  the  choice  of  laser.  The  choice  of  fluorescence  wavelength  is  driven  by  the  choice 
of  detector.  The  radiative  lifetime  of  the  upper  state  limits  the  vertical  resolution  of  the 
lidar.  Nonradiative  quenching  of  the  upper  state  decreases  the  number  of  photons  that 
would  otherwise  be  produced  by  fluorescence. 

The  upper  rotational  level  and  the  lower  rotational  level  of  the  fluorescence  transition 
are  almost  determined  after  the  lower  rotational  level  of  the  excitation  transtion  is  chosen. 
The  rotational  transition  must  obey  the  selection  rule  A  J  =  0,±1,  with  the  condition 
that  A  J  =  0  is  forbidden  for  electronic  systems  in  which  both  states  have  zero  angular 
momentum  along  the  internuclear  axis.1  The  transition  J  =  0  — ►  J  =  0  is  always 
forbidden.  Transitions  with  A  J  =  Jupp4r  -  Ji0VltT  =  -1  are  in  the  P  branch,  A  J  =  +1  are 
in  the  R  branch,  and  A  J  —  0  are  in  the  Q  branch. 

After  the  upper  rotational  level  of  the  excitation  transition  is  chosen  from  two,  or  in 
some  cases  three,  possible  choices,  the  lower  rotational  level  of  the  fluorescence  transition 
can  be  chosen  so  that  the  lidar  detects  any  combination  of  the  two  (or  three)  fluorescing 
branches.  The  filter  in  the  lidar  receiver  determines  which  brancnes  are  detected.  If  the 
bandwidth  of  the  filter  is  wide  enough  (perhaps  on  the  order  of  10-30  A)  then  all  branches 
can  be  detected  and  the  signal  is  two  (or  three)  times  higher  than  it  would  otherwise  be 
if  detecting  only  one  branch.  However,  there  is  a  tradeoff  because  a  wide  bandwidth  filter 
allows  more  background  radiation  to  be  detected. 

Figure  3  shows  the  four  lidar  scenarios  along  with  the  list  of  lasers  and  the  specific 
electron/vibrational/rotational  transitions  that  we  consider  for  each.  The  excitation  and 
fluorescence  wavelengths  are  also  shown.  The  wavelengths  are  calculated  using  spectro¬ 
scopic  constants  from  Ref.  [7].  Appendix  Appendix  A  describes  these  calculations  in  detail. 

1In  the  parlance  of  molecular  spectroscopy  such  a  state  is  designated  by  E.  For  example,  the  1**  negative 
system  refers  to  transtions  <Tom  Nj  (.X^E+)  to  Nj(B2E+).  Both  states  are  E  states  and  therefore  A  J  0. 
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Jbkgrnd  II. 

1)  LIF 

4)  Pump/LIF 

N2(B) 


1 


Ns(A)J 


1st  pos 


1)  Ti-Saph.dye:  v-O,J-0  — >  v-2,J-5  — >  v-O.J-4,5,0 

(7730.5)  (7717.SU 


N2:  v-0  — >  v-2  — >  v-0 
(7737)  (7737) 


2)  N2:  v-0  —  >  v-0  — >  v-0 
(3371)  (3371) 

dye:  v-O.J-0  — >  v-O.J-5  —  >  v-O.J-4,5.0 
(3371.9)  (33C79) 

(3309.7) 

(3371.9) 


Nd:YAG:  V-5.J-12  —  >  V-4.J-13  -->  v-2.J-12.13, 14 
(10047)  (7445.9) 

(74*®;|{ 


3)  dye:  v-I.J-0  —  >  v-O.J-5  — >  v-O.J-4,0 
(3914.3)  (3907.8) 

(3914.3) 


(7488. 


4)  pump:  dye:  v-O.J-0  — >  v-O.J-5 

(10470.3) 

Nd:YLF:  v-O.J-3  —  >  v-O.J-2 
(10472.0) 

LIF:  N2:  v-0  — >  v-0  — >  v-0 
(3371)  (3371) 


Figure  3:  Specific  transitions  and  the  corresponding  wavelengths  for  the  four  lidar  scenar¬ 
ios. 
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2.2  LIF  Lidar  Signal 


For  the  conventional  LIF  lidar  the  number  of  photons  n*,(r)  detected  per  laser  pulse  by 
the  lidar  receiver  in  a  time  period  At,  at  a  time  t  =  2 r/c  after  the  start  of  a  laser  pulse  is 


p  \  £  n2  av  ~v"o 

n„,(r)  =  vK^K^Nt'VRt'V(J)GF:::;„'  W Tu{r)Td(r)ArT{r). 


This  equation  assumes:  1)  the  fluorescence  medium  is  not  saturated,  and  2)  the  spec¬ 
tral  shape  of  the  laser  radiation  is  approximately  the  same  as  the  spectral  shape  of  the 
excitation. 


Each  factor  in  Eq.  2  is  described  below. 


•  q  is  the  quantum  efficiency  of  the  detector. 

•  K  is  the  transmissivity  of  the  receiver  optics. 

•  Ea  is  the  laser  energy  per  pulse. 

•  Dr  is  the  diameter  of  the  telescope  aperture. 

•  Nti„  is  the  density  of  N2  or  N2  in  electronic  state  e  and  vibrational  state  v  (the 
species  to  be  excited). 

•  Re  v(J)  is  the  fraction  of  the  Ne>v  molecules  which  are  in  rotational  state  J  (see  Eq.  1 
above). 

•  Fe'-e"  IS  the  Franck-Condon  factor  associated  with  a  transition  from  electronic/ vibrational 
state  e'/v'  to  electronic /vibrational  state  e"/v"  (i.e.,  the  probability  that  N2(e',t/) 
decays  to  N2(e",  v")  instead  of  other  vibrational  states). 

•  G  is  the  fraction  of  fluorescence  radiation  that  is  within  the  bandwidth  of  the  receiver. 


•  Og!Zf>>  is  the  absorption  cross  section  for  absorption  of  a  photon  by  a  molecule  in  the 
electronic/vibrational  state  en/v"  to  produce  a  molecule  in  the  electronic/vibrational 
state  e'/v'.  The  absorption  cross  section  is  calculated  in  the  following  way: 


c =  — — 
BnAu 


(mks  units), 


where 


-  A  is  the  transition  probability, 
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-  A  is  the  wavelength  of  the  transition,  and 

-  A  u  is  the  bandwidth  of  the  transition  (primarily  doppler). 

•  Q  is  the  quenching  factor  for  nonradiative  decay  of  the  excited  molecule. 

•  The  factor  4jt  indicates  that  fluorescence  emission  is  assumed  isotropic. 

•  Tu(r)  and  Td(r)  are  the  transmissivities  of  the  atmosphere  from  the  lidar  up  to 
range  r  at  the  wavelengths  of  the  upwardly  propagating  and  downwardly  propagating 
radiation,  respectively. 

•  A r  =  c  At/2,  where  At  is  the  integration  time  of  the  detector. 

•  r(r)  is  the  fluorescence  lifetime  correction  factor  [6).  This  factor,  which  is  a  dimen¬ 
sionless  quantity  less  than  or  equal  to  unity,  is  present  due  to  the  nonzero  lifetime  of 
the  excited  state  and  the  fact  that  the  lidar  return  at  time  t  =  2 r/c  may  be  due  to 
photons  from  ranges  less  than  r  as  well.  This  factor  approaches  nr  :ty  as  the  distance 
into  the  fluorescence  medium  increases.  For  fluorescence  lifetimes  short  compared  to 
detector  integration  time  this  factor  is  near  unity  for  all  range  points. 

Section  5  presents  lidar  parameters  and  results  of  calculations  for  numbers  of  photons 
detected  by  hypothetical  LIF  lidars  (see  Table  6). 

3  Pump/LIF  Lidar 

3.1  General  Description 

The  pump/LIF  lidar  is  similar  to  a  conventional  LIF  lidar.  However  the  excited  species 
that  it  probes,  in  an  LIF  lidar  manner,  is  created  by  first  laser  pumping  from  a  lower  state. 
This  lidar  utilizes  two  lasers.  The  first  laser,  the  so-called  pump  laser,  excites  nitrogen  from 
Nj(j43£J)  to  N2(BsII,)  (1**  positive  system).  The  second  laser,  the  so-called  the  probe 
laser,  excites  nitrogen  from  Nj(.Bsn,)  to  Nj(C3IIu)  (2nd  positive  system).  The  receiver 
detects  fluorescence  from  Nj(C3IIu)  to  the  N2(JB3rij). 

The  impetus  for  this  type  of  lidar  is  the  desirability  of  doing  LIF  lidar  on  the  2nd 
positive  system.  Such  a  lidar  detects  radiation  in  the  near  UV  for  which  higher  quantum 
efficiency  photomultiplier  tubes  exist.  L  ^fortunately,  the  density  of  Nj(BsII,)  that  exists 
during  moderate  aurora  is  too  low  for  this  type  of  lidar  by  itself  to  be  feasible.  However, 
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the  density  of  Nj(A*E+)  is  quite  high.  Therefore,  another  laser  is  used  to  create  Nj(BsIIf) 
from  the  abundance  of  Nj(AsE+). 

Another  advantage  of  the  pump/LIF  lidar  becomes  evident  after  calculating  the  number 
of  photons  expected  to  be  detected.  The  number  of  photons  detected  is  proportional  to 
the  inverse  of  the  square  of  the  laser  divergence,  if  the  divergences  of  both  lasers  are  the 
same.  (If  the  divergences  of  the  pump  and  probe  lasers  are  different,  then  take  the  larger 
of  the  two  divergences).2  This  property  of  the  pump/LIF  lidar  allows  for  the  increase  of 
the  signal  by  placing  beam  expanders  on  the  outputs  of  the  lasers.  The  only  limitation  is 
the  diffraction  limit  of  the  beam  expanders  (i.e.,  the  size  of  the  beam  expander  aper 

Similar  issues  guide  the  choice  of  transitions  for  the  pump/LIF  lidar  as  those  that  gu.ae 
the  choice  for  the  conventional  LIF  lidar.  The  pump/LIF  lidar  has  the  pump  transition 
(induced  by  the  pump  laser),  the  excitation  transition  (induced  by  the  probe  laser),  and  the 
fluorescence  transition  (detected  by  the  receiver).  Figure  3  shows  the  specific  transitions 
that  we  consider.  In  that  figure  the  pump/LIF  lidar  is  depicted  as  scenario  number  four. 


3.2  Pump/LIF  Lidar  Signal 


For  the  pump/LIF  lidar  the  number  of  photons  »„-,(r)  detected  per  laser  pulse  by  the  lidar 
receiver  in  a  time  period  At,  at  a  time  t  =  2 r/c  after  the  start  of  a  laser  pulse  is 


n*t(r)  =  VK- 


!  W 


^~^^NA,.y^-’’Fg^"G^^^-Tr(r)T,(r)Ur)T(r)Ar 

(4) 

In  Eq.  4  A,  B,  and  C  refer  to  the  three  electronic  levels  involved  in  the  pump/LIF  process. 
The  pump  laser  pumps  from  A  to  B.  The  probe  laser  excites  from  B  to  C.  The  receiver 
detects  fluorescence  radation  from  the  decay  of  C  to  B.  Each  factor  in  Eq.  4  which  does 
not  also  occur  in  Eq.  2  is  described  below.  Factors  which  are  also  in  Eq.  2  are  described 
in  the  text  after  that  equation. 


•  0i  is  the  laser  divergence.  Equation  4  assumes  both  lasers  to  matched  (i.e.,  same 
divergences). 

•  ea  is  the  laser  energy  per  pulse  of  the  pump  laser. 

•  Eb  is  the  laser  energy  per  pulse  of  the  probe  laser. 

3For  a  conventional  LIF  lidar  the  number  of  photons  detected  is  independent  of  the  laser  divergence  as 
long  as  the  receiver  field  of  view  completely  encompasses  the  laser  beam. 
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•  A ab  is  the  wavelength  corresponding  to  the  pump  transition. 

•  A xb  is  the  wavelength  corresponding  to  the  fluorescence  transition. 

•  Na.v,  is  the  density  of  Nj(A,  vp),  the  state  excited  by  the  pump  laser.  vp  is  the 
vibrational  level. 

•  Ra,v,(Jp)  is  the  fraction  of  the  Na,v,  molecules  which  are  in  rotational  state  JP  (see 
Eq.  1  above). 

•  is  the  absorption  cross  section  for  the  pump  transition.  The  absorption  cross 
section  is  calculated  using  Eq.  3. 

•  Obcv"  is  the  absorption  cross  section  for  the  excitation  transition. 

•  Fbcv"  is  the  Franck-Condon  factor  associated  with  fluorescence  transition  vibrational 
band. 

•  G  is  the  fraction  of  fluorescence  radiation  that  is  within  the  bandwidth  of  the  receiver. 

•  Tp(r),  Ti  (r) ,  and  Tu(r)  are  the  transmissivities  of  the  atmosphere  from  the  lidar  up 
to  range  r  for  the  pump,  excitation,  and  fluorescence  wavelengths,  respectively. 

Section  5  presents  lidar  parameters  and  results  of  calculations  for  numbers  of  photons 
detected  by  hypothetical  pump/LIF  lidars  (see  Table  6). 

4  Density  of  Excited  and  Ionized  Molecular  Nitrogen 

The  calculation  of  LIF  or  pump/LIF  lidar  signals  requires  knowledge  of  the  density  of  the 
specific  vibrational  state  that  is  to  be  excited  (LIF  lidar)  or  pumped  (pump/LIF  lidar). 
The  fraction  of  these  that  are  in  a  specific  rotational  state  is  then  calculated  with  Eq.  1. 
The  density  of  a  vibrational  state  that  radiatively  decays  to  other  states  can  be  derived 
from  observed  column  emission  rates,  of  the  relevant  bands,  from  aurora.  The  density  of 
a  state  that  does  not  radiatively  decay  (e.g.,  (A'JE^))  requires  a  different  technique, 

which  we  describe  below.  Of  all  the  lidar  scenarios  that  we  consider  (see  Figure  3)  only 
one,  the  1**  negative  LIF  lidar,  requires  knowledge  of  the  density  of  a  nonradiative  state. 
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4.1  Density  of  Excited  N2 


The  1**  positive,  2nd  positive,  and  pump/LIF  lidars  depicted  in  Figure  3  involve  exciting 
Nj(AsE+)  or  Nj(B*II#).  The  densities  of  specific  vibrational  states  of  Nj(AsE+)  can  be 
derived  from  observed  column  emission  rates  of  the  relevant  bands  in  the  Vegard-Kaplan 
system.  The  densities  of  specific  vibrational  states  of  Nj(5sIl#)  can  be  derived  from 
observed  column  emission  rates  of  the  relevant  bands  in  the  1*‘  postive  system.  Column 
emission  rates  for  type  IBC3  aurora  are  reported  in  Reference  [5].  The  relevant  bands 
are  all  those  which  have  an  upper  vibrational  level  identical  to  the  vibrational  level  of  the 
molecule  to  be  excited. 

For  steady  state  conditions  the  density  Ne>v  of  Nj  in  electronic  state  e  and  vibrational 
state  v  is 


where 

•  is  the  column  emission  rate  (in  Rayleighs)  due  to  the  radiative  decay  of  the 
electronic/ vibration  state  ejv  to  the  electronic/vibrational  state  e'/v', 

•  rC|U  is  the  radiative  lifetime  of  the  tfv  state,  and 

•  T  is  the  thickness  of  the  aurora,  which  we  take  as  10  km. 

The  quantity  in  parentheses  is  the  volume  emission  rate.  Note  that  this  calculation  does 
not  require  knowledge  of  lifetimes  associated  with  non-radiative  decay  mechanisms,  which 
are  typically  difficult  to  estimate. 

Table  1  shows  the  density  of  the  Nj  states  needed  to  calculate  numbers  of  photons 
for  the  lidar  scenarios  that  we  consider.  Also  shown  are  the  transitions  that  are  used 
to  determine  the  densities,  the  corresponding  observed  column  emission  rates,  and  the 
radiative  lifetimes  of  the  states. 

4.2  Density  of  N£ 

Here  we  estimate  the  density  of  N,  (X3E^,v=0)  since  this  is  the  specific  vibrational  state 
to  be  excited  with  the  l*1  negative  LIF  lidar  that  is  depicted  in  Figure  3. 

Since  Nj  (XaE^",v=0)  does  not  radiatively  decay  we  cannot  use  the  above  technique  to 
calculate  its  density.  Instead,  we  estimate  the  ratio  of  the  density  of  N,  (X2E+,v=0)  to  the 
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LIF 

excitation 

transition1 

relevant 

auroral 

transition1 

column 
emission 
rate  (kR) 

radiative 

lifetime 

(sec) 

density  (cm-3) 
(lower  state 
of  LIF  trans.)1 

VK:  (A,v  =  0)  ->  (X,  all  v) 

w 

2 

111 

P  f  f  Jus 

VK:  (A,  v  =  5)  — ► (X,  all  v) 

mm 

2 

IP:  (B,v  =  0)  -+  (A,  all  v) 

9xl0~® 

■‘ill  states  are  of  Nj  2 assume  10  km  aurora  thickness 


Table  1:  Densities  of  Nj  vibrational  states  to  be  excited  by  the  considered  lidars.  All 
parameters  used  to  calculate  the  densities  are  included. 

density  of  Nj  (J93E+)  and  then  calculate  the  density  of  Nj  (£*£+),  which  does  radiatively 
decay. 

In  the  aurora  Nj  is  produced  primarily  by  electron  impact  ionization  of  ground  state 
Nj.  Nj  losses  are  due  primarily  to  reactions  with  O  and  Oj  (see  Table  2).  In  steady  state 
the  density  Nx+  of -Nj(X*E^)  is  then  simply  product  of  the  electron  impact  ionization 
production  rate  and  the  loss  time  due  to  all  mechanisms:3 

JV,+  =  N0{a$v)  r,+  (6) 

where 

•  N0  is  the  density  of  ground  state  Nj, 

•  N0(al+v)  is  the  production  rate  due  to  electron  impact  ionization  of  ground  state  Nj 
(v  is  the  relative  speed  between  the  impacting  electron  and  the  molecule,  a*+  is  the 
differential  cross  section,  a  function  of  v,  and  (•  •  •)  denotes  an  average  with  respect 
to  the  velocity  destribution  function),  and 

•  rz+  is  the  loss  time  of  Nj  (X*E* )  due  to  all  loss  mechanisms. 

The  density  Nz+>v=0  of  the  specific  vibrational  state  Nj  (X*E+,v=0)  is  then  approximately 

N,+'V=0  =  N,+F?;°_s  (7) 

where  is  the  Franck-Condon  factor  for  the  transition  NJ(X1E^,v=0) — ►  Nj (X*E+,v=0).4 

3This  disregards  spatial  tranport  processes  and  production  due  to  electron  impact  of  other  states  besides 
ground  state  N2.  We  assume  these  two  effects  to  be  small. 

disregards  production  of  Nj uue  to  electron  impact  with  N,  (X,Ej,v#0). 
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The  density  JV»+  of  Nj  (B*EJ)  is  approximately 


JV»+  =  N0(al+v)rh+ 


(8) 


where  each  variable  is  analogous  to  the  corresponding  one  in  Eq.  6.  r»+ ,  the  loss  time  of 
Nj  (£*££),  is  due  primarily  to  radiative  decay. 

Since  (£*££)  radiatively  decays  its  density  is  also  determined  by  the  technique 
described  in  the  last  section.  Therefore,  N^+  is  also  written  as 


(») 


where  the  emission  is  considered  to  be  from  the  entire  system. 
Dividing  Eq.  6  by  Eq.  8  and  substituting  in  Eq.  9  we  get 


^*+,«=o  “ 


T' 


<T>  ■*'*+-* 


(10) 


The  following  values  are  used  for  the  quantities  on  the  right  hand  side  of  Eq.  10: 


(<rx+v)  /  (^b+«) 
r,+ 


rB+ 

Eg+ 

T 

jpO-O 
*  *+-* 


4.3  (calculated  using  diff.  cross  sections  from  Ref.  [8]) 

4.4  msec  (at  100  km) (see  Table  2  below) 

90  nsec  (primarily  radiative) 

80  kiloRay’eighs 

10  km 

0.09 


The  density  of  N,  (XaE+,v=0)  is  then 


#*+,v=o  =  1.4  x  10*  cm"*. 


5  Signal,  Noise,  and  Signal  to  Noise  Ratio 

5.1  Signal 

Here  we  present  the  results  of  calculations  for  the  number  of  signal  photons  from  100  km 
for  the  four  lidar  scenarios  depicted  in  Figure  3.  The  numbers  of  photons  for  the  LIF  lidars 
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reaction 

reaction  time 

1.  Nf  +  O,  — *  W,  +  0+ 

2.  Nf+0—>NO++N 

3.  W,+  +  O  — »  Nt  +  O* 

4.  Nf  +  e~  — *  N  +  N 

6.50  msec 
14.5  msec 
212  msec 

900  sec 

Table  2:  Reactions  primarily  responsible  for  decay  of  Nf  at  100  km  (T=200°K),  and  their 
reaction  times.  (From  Ref.  9.) 

are  calculated  using  Eq.  2.  The  number  of  photons  for  the  pump/LIF  lidar  is  calculated 
using  Eq.  4.  The  lidar  parameters  considered  for  these  calculations  are  listed  in  Table  3. 

The  relevant  absorption  cross  sections  required  for  these  calculations  are  listed  in  Tar 
ble  4.  The  absorption  cross  sections  are  calculated  using  Eq.  3.  Also  shown  in  the  table 
are  transitions  probabilities  and  linewidths.  The  linewidths  are  assumed  to  be  due  pri¬ 
marily  to  Doppler  broadening.  The  temperature  at  100  km  is  taken  to  be  200°K  [11].  The 
transition  probabilities  are  obtained  from  Ref.  [12]. 

The  relevant  Franck-Condon  factors  required  for  these  calculations  are  listed  in  Table  5. 
These  values  are  obtained  from  Ref.  [12]. 

Table  6  shows  the  number  of  photons  detected  per  laser  pulse  in  a  10  /isec  time  period 
from  100  km  altitude  during  an  IBC3  aurora  for  each  of  the  lidar  scenarios  that  we  consider. 

5.2  Uncertainties  in  the  Signal 

We  consider  four  noise  sources:  1)  signal  shot  noise,  2)  detector  noise,  3)  sky  background, 
and  4)  auroral  background  emission.  The  uncertainty  (standard  deviation)  in  the  number 
of  signal  photons  nMf(r)  per  laser  pulse  in  one  detector  integration  time  period  At  is 

=  ni//(r)- 

The  uncertainty  in  the  number  of  effective  photons  in  a  time  period  At  due  to 
detector  noise  is  o<ut  —  nxj£.  For  a  photomultiplier  tube  n<ut  can  be  calculated  from  the 
dark  current  tj,  which  is  a  quantity  typically  specified  with  a  tube: 

»*.  =  ^  (PMT)  (II) 

where  e  is  the  electron  charge  and  G  is  the  gain  of  the  dynode  chain. 

The  uncertainty  in  the  number  of  sky  background  photons  is  o(*v  =  n)£.  The 
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E0  = 

100  mJ  (dye  laser) 

10  mJ  (nitrogen  laser) 

2  J  (Nd:YAG  laser) 

2.5  J  (Ti:Sapphire  laser) 

1J  (NdrYLF  laser) 

Dr  — 

1  m 

9r  = 

0.5  mrad 

At  = 

10  /isec 

AA  = 

10  A 

V  = 

0.3  (1*‘  neg,  2nd  pos) 

0.1  (1*  pos) 

K  = 

0.2  (1**  neg,  2nd  pos) 

0.3  (l*1  pos) 

r  = 

100  km 

T(r)  = 

0.28  (2nd  pos)t 

0.5  (1**  neg)* 

0.8  (1**  pos)f 

*  Calculated  with  LOWTRAN  [10] 


Table  3:  Lidar  parameters  used  for  the  lidar  signal  calculations. 

number  of  sky  background  photons  detected  in  a  time  period  At  is 

«,*  =  riK^AHAXAt  (12) 

he 

where  B\  is  the  sky  background  spectral  radiance  at  wavelength  A,  ft  is  the  receiver  field 
of  view  solid  angle,  A  A  is  the  receiver  bandwidth,  r}  is  the  detector  quantum  efficiency,  and 
K  is  the  transmissivity  of  the  receiver  optics.  The  lidars  we  consider  can,  at  best,  operate 
during  a  moonlit  night  for  which  B  =  2.5  x  10-9  W/m*-ster-A  [13]. 

Nitrogen  molecules  which  are  excited  directly,  by  auroral  electron  impact,  to  the  upper 
state  of  the  fluorescence  transition  may  radiatively  decay  and  be  detected  by  the  lidar 


1*‘  pos  (2-0) 

iUffgltMl 

sgisaiiiiii 

l*1  neg  (1-0) 

A  (sec-1) 
Ai/(200°K)  (GHz) 
oa  (cm*) 

4.4  xlO4 
0.74 

1.43X10-14 

7.84  x10s 
0.54 

6.55  XlO-18 

1.39X107 

1.7 

3.70x10" 18 

9.64x10® 

1.5 

3.92  XlO-18 

Table  4:  Transition  probabilities,  doppler  linewidths,  and  photoabsorption  cross  sections 
used  for  the  lidar  signal  calculations. 
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1**  pos  (2-0) 

1**  pos  (4-2) 

2nd  pos  (0-0) 

1*‘  neg  (0-0) 

0.1617 

0.2976 

0.4527 

0.6481 

Table  5:  Franck-Condon  factors  used  for  lidar  signal  calculations. 


laser 

LIF 

1**  pos 

LIF 

2nd  pos  (0-0) 

LIF 

1**  neg  (0-0) 

pump  (1**  pos  (0-0))- 
LIF  (2nd  pos  (0-0)) 

nitrogen: 

0.019 

9x10"* 

0.0024 

- 

dye: 

1.7xl0~6 

0.0047 

0.30* 

Nd:YAG: 

- 

- 

- 

Nd:YLF: 

- 

- 

- 

3.0* 

Ti:Sapphire: 

1.4 

- 

0.0094 

- 

^Numbers  appear  in  pump  laser  row.  LIF  laser  is  nitrogen. 


Table  6:  Number  of  photons  detected  per  laser  pulse  for  all  lidars  considered. 

receiver.  These  photons  constitute  unwanted  background,  which  we  refer  to  as  auroral 
background  emission.  The  uncertainty  in  the  number  of  auroral  background  emission 
photons  naur  is  oaur  =  nlJJr.  The  number  of  auroral  background  emission  photons  detected 
in  a  time  period  At  is 

>W  =  vK^E’;:;:a^DlAtTd{r)  (13) 

4  ID 

where  is  the  column  emission  rate  (in  Rayleighs)  due  to  transitions  from  the  elec¬ 

tronic/vibrational  state  e'/v'  to  cn/vn,  G  is  the  fraction  of  these  photons  that  is  within  the 
bandwidth  of  the  receiver  (see  discussion  in  Section  2),  and  all  other  quantities  have  the 
usual  definitions. 

Table  7  below  shows  the  the  number  of  sky  background  photons  (moonlit  night)  and 
the  number  of  auroral  background  emission  photons  detected  in  a  time  At. 


l*f  pos  (2-0) 

1**  pos  (4-2) 

EjESCgg 

1**  neg  (0-0) 

auroral  emission*  (kR) 

53.8 

24.4 

50 

auroral  emission  (#  photons) 

■ 

0.31 

0.74 

sky  bkgrnd  (#  photons) 

0.009 

0.008 

Table  7:  Number  of  sky  background  photons  and  auroral  background  emission  photons  in 
time  At. 
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If  all  noise  processes  are  statistically  independent  of  each  other  then  the  total  uncer¬ 
tainty  on  in  the  number  of  photons  for  one  laser  pulse  is 

On  =  +  °)ky  +  al,t  +  o\w  =  +  ntky  +  n^t  +  naur  (14) 

For  m  laser  pulses  each  photon  number  is  multiplied  m. 


5.3  Signal/Noise 

The  signal  to  noise  ratio  ( S/N)m  for  m  laser  pulses  is 


(£) 

\NJm  y/nHl{r) 


(15) 


r)  +  +  n<ut  +  flaw 

The  following  table  shows  the  signal  to  noise  ratios  (for  one  laser  pulse),  for  all  the  lidars 
considered,  using  Eq.  15  and  the  information  and  parameters  presented  in  this  section. 


laser 

LIF 

1*‘  pos 

LIF 

2nd  pos  (0-0) 

LIF 

1*‘  neg  (0-0) 

pump  (1**  pos  (0-0))- 
LIF  ( 2nd  pos  (0-0)) 

nitrogen: 

2xl0'6 

- 

dye: 

0.087 

3.0xl0“6 

1 

0.38* 

Nd:YAG: 

HUB 

- 

- 

- 

Nd:YLF: 

- 

- 

- 

1.6* 

Ti:Sapphire: 

i.i 

- 

0.011 

- 

*  Numbers  appear  in  pump  laser  row.  LIF  laser  is  nitrogen. 


Table  8:  Signal  to  noise  ratios  for  one  laser  pulse  (S/N) *  for  all  the  lidars  considered. 


The  signal  to  noise  ratio  scales  with  the  parameters  in  the  following  way: 

•  Auroral  background  fluorescence  limited: 

(D^oc  (m'lA-Ar)1'IA£0£>r 

•  Shot  noise  limited: 

(f)m  oc  (mriK  ArXEoNt'u)1^2  Dr 
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6  Conclusions 


We  have  described  two  types  of  molecular  nitrogen  fluorescence  lidars  which  can  be  used 
to  determine  exciteu  or  ionized  Nj  density  in  the  auroral  ionsphere.  In  addition,  the  lidars 
can  be  used  for  high  spatial  resolution  diagnostics  during  ionospheric  high  power  radio 
frequency  heating  experiments  which  are  often  performed  during  moderate  aurora.  The 
first  lidar  is  a  conventional  LIF  lidar  which  operates  in  the  first  postive  system  of  N2  or 
the  first  negative  system  of  Nj.  Several  different  lasers  are  considered.  The  choice  of  laser 
determines  the  bands  of  excitation  and  fluorescence.  The  second  lidar  is  a  pump/LIF  lidar. 
Two  lasers  are  utilized.  The  first  laser  pumps  N2  from  A3£+  to  BSII#  ((0-0)  band,  10477 
A).  The  second  laser  then  does  LIF  on  the  second  positive  system  ((0-0)  band,  3371  A). 

We  have  calculated  the  signal  to  noise  for  both  types  of  lidars  using  reasonable  lidar 
parameters.  We  have  considered  four  noise  processes:  1)  signal  shot  noise,  2)  auroral  back¬ 
ground  emission,  3)  ambient  sky  background,  and  4)  detector  noise.  We  have  considered 
several  different  electronic  systems  of  N2  and  Nf  (l4t  pos,  2nd  pos,  1**  neg)  and  several 
different  lasers. 

The  most  promising  lidars  are 

1.  A  LIF  lidar  operating  in  the  l4*  positive  system  (2-0)  band,  using  a  tunable  Ti- 
Sapphire  laser.  The  signal  to  noise  for  one  laser  pulse  is  [S/N)\  =  1.1.  Required 
averaging  time  for  S/N  =  100  is  165  sec  (at  50  Hz). 

2.  A  LIF  lidar  operation  in  the  l4t  positive  system  (4-5)  band,  using  a  seeded  Nd:YAG 
laser  (tunable  over  «  1  A).  For  this  lidar  (S/N)i  =  0.24.  Required  averaging  time 
for  S/N  =  100  is  5790  sec  (at  30  Hz). 

3.  A  pump/LIF  lidar  using  a  Nd:YLF  laser  as  a  pump  in  the  14‘  positive  system  and  a 
dye  laser  for  LIF  in  the  2nd  positive  system.  For  this  lidar  {S/N)  i  =  1.6.  Required 
averaging  time  for  S/N  —  100  is  130  sec  (at  30  Hz). 

The  pump/LIF  lidar  is  preferable  to  the  conventional  LIF  lidar  since  higher  signal  to 
noise  ratios  can  be  achieved.  The  primary  reason  for  this  is  that  the  signal  is  proportional 
to  the  square  of  the  inverse  of  the  divergences  of  the  lasers.  The  divergences  can  be  made 
very  small  with  the  use  of  beam  expanders. 

The  conventional  LIF  lidar  is  preferrable  in  the  short  term  since  it  utilizes  only  a  single 
laser.  It  is  easier  to  adapt  an  existing  Rayleigh  lidar  to  become  a  LIF  lidar.  Utilizing  a  dye 
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laser  for  this  type  of  lidar  is  preferable  since  the  highest  signal  to  noise  ratio  can  be  achieved. 
However,  since  many  Rayleigh  lidars  use  high  power  Nd:YAG  lasers,  it  is  straightforward 
to  convert  such  a  lidar  into  a  LIF  lidar  as  described  in  this  report.  All  that  is  required  is  a 
new  filter  and  a  new  photomulitplier  tube  in  the  receiver.  The  nitrogen  laser  approach  for 
the  conventional  LIF  lidar  is  attractive  since  no  tuning  is  required.  However,  this  desirable 
feature  is  offset  by  the  low  power.  There  is  a  question  as  to  availability  of  a  first  positive 
nitrogen  laser. 
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Appendix  A  Wavelength  Calculations 


All  wavelengths  presented  in  this  report  axe  calculated  theoretically  in  the  manner  de¬ 
scribed  here. 

A  rotational  level  of  a  diatomic  molecule  such  as  JVj  and  Nf  is  specified  by  its  term 
energy  T  =  E/hc,  where  E  is  the  energy.  The  term  energy  of  a  rotational  level  can  be 
written  as 

T  =  T0  +  G{v)  +  FV(J)  (16) 

where  T0  is  the  electronic  term,  G(v)  is  the  vibrational  term,  and  FU(J)  rotational  term, 
v  and  J  are  the  vibrational  and  rotational  quantum  numbers,  respectively.  In  general,  the 
rotational  term  is  a  function  of  the  vibrational  quantum  number  as  well  as  the  rotational 
quantum  number. 

The  vibrational  and  rotational  terms  are  written  as  series  in  the  vibrational  and  rota¬ 
tional  quantum  numbers  respectively: 


G(v)  =  u>e(v  +  1/2)  -  uexe(v  +  1/2)2  -I-  weye(v  +  l/2)s  4-  u>eze(v  +  1/2)4  -f - (17) 

FV(J)  =  BvJ(J  +  l)~DvJ*(J  +  l)a  +  HvJa(J  +  l)a  +  ---  (18) 

where 

Bv  =  Be  -  cte{v  +  1/2)  +  7e(v  +  1/2)J  +  6t(v  +  l/2)s  +  •  •  •  (19) 

Dv  —  Dt  +  /3«(v  +  1/2)  +  •  •  •  (20) 

Hv  =  He+  -  (21) 

The  quantities 

^et  %et  !/<>  ^ei  D e,  Bet  ®«»  7«)  ^tt  fit 


are  spectroscopic  constants  which  can  be  looked  up  in  tables.  We  use  tables  from  Ref.  [7]. 

The  wavelength  of  radiation  resulting  from  a  radiative  decay  between  two  rotational 
levels  is 

Tu-Ti 

where  Tu  and  Tj  are  the  term  energies  for  the  upper  and  lower  transition  levels  respectively. 
All  wavelengths  in  this  report  are  calculated  using  Eq.  22  with  Eqs.  16  to  18.  Calculated 
wavelengths  have  uncertainties  due  to  uncertaintes  in  spectroscopic  constants  and  unavail¬ 
ability  of  some  constants.  Uncertainties  may  be,  at  most,  on  the  order  of  a  few  tenths  of 
an  Angstrom. 
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Appendix  B  Tests  at  UCLA  Plasma  Physics  Labora¬ 
tory 

B.l  Introduction 

In  February,  1993  one  of  the  authors  (Gamer)  spent  two  weeks  at  the  UCLA  Plasma 
Physics  Laboratory  (UCLA  PPL)  as  part  of  a  collaborative  effort  between  Phillips  Lab¬ 
oratory  (Geophysics  Directorate,  GPIM  branch)  and  UCLA  PPL  to  study  the  feasibility 
of  a  molecular  nitrogen  fluorescence  lidar  as  a  diagnostic  for  high  power  rf  heating  experi¬ 
ments  of  the  ionosphere.  The  specific  purpose  of  the  visit  was  to  investigate  laser  induced 
fluorescence  of  the  first  positive  system  of  N3  in  a  laboratory  based  nitrogen  plasma  to 
determine  its  feasibility  as  the  basis  for  a  lidar.  The  GPIM  branch  supplied  a  dye  laser 
for  these  studies.  UCLA  PPL  personnel  directly  involved  in  the  effort  were  Dr.  Ralph 
Wuerker,  Mr.  Tetsuo  Fukuchi,  and  Dr.  Alfred  Wong.  The  general  laboratory  facilities  and 
support  personnel  were  available  as  well. 

The  results  of  the  experiments  were  inconclusive.  We  were  not  successful  in  demon¬ 
strating  laser  induced  fluorescence  of  the  first  positive  system.  We  discuss  possible  reasons 
for  this  below.  Additional  experimentation  is  required. 

The  UCLA  PPL  has  an  ongoing  research  program  aimed  at  investigating  the  spec¬ 
troscopy  of  Nj  and  N,  in  laboratory  simulations  of  ionospheric  plasmas.  During  this 
research  they  have  successfully  demonstrated  laser  induced  fluorescence  in  the  the  second 
positive  system  of  Nj  and  the  first  negative  system  of  N3  [14].  They  do  not  have  a  laser 
suitable  for  excitation  in  the  first  positive  system.  For  the  second  positive  system  they 
used  a  nitrogen  laser  operating  in  the  (0-0)  band  at  3371  A  and  detected  fluorescence  in 
the  same  band  and  from  several  other  bands  with  upper  vibrational  level  v  =  0.  For  the 
first  negative  system  they  used,  in  different  experiments,  a  nitrogen  laser  and  a  dye  laser 
laser.  Both  lasers  operated  in  the  in  the  (1-0)  band  at  4278  A.  The  dye  laser  operated 
at  a  specific  rotational  level  within  the  band.  In  these  LIF  experiments  they  measured 
lifetimes  of  the  observed  fluorescence  transitions  and  absorption  cross  sections  of  the  exci¬ 
tation  transitions.  For  the  first  negative  system  they  compared  the  relative  merits  of  wide 
band  excitation  (nitrogen  laser)  and  narrow  band  excitation  (dye  laser). 

Our  calculations,  presented  in  this  report,  indicate  that  the  first  positive  system  of  Nj 
may  be  more  promising  for  lidar  use  than  the  other  systems.  Therefore,  we  set  out,  with  the 
UCLA  PPL  group,  to  demonstrate  first  postive  laser  induced  fluorescence.  Specifically,  we 
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set  out  to  demonstrate  the  first  scenario  depicted  in  Figure  3.  In  one  part  of  that  scenario  a 
dye  laser  is  used  to  excite  the  transition  Nj(AsE+,v=0,J=6) — ♦  Nj(Bsn#,v=2,J=5).  The 
wavelength  is  theoretically  (see  Appendix  Appendix  A)  7737  A.  Fluorescence  radiation  is 
then  detected  from  one  or  all  three  of  the  branches  of  the  reverse  transition.  The  (2-0) 
band  of  the  first  positive  system  offers  the  best  combination  of  wavelength,  cross  section, 
and  Franck-Condon  factor  for  lidar  use  (see  discussion  in  Section  2.1).  The  J  =  5  rotational 
level  is  expected  to  be  the  most  abundant  rotational  level  at  100  km  (200° K). 

The  specific  goals  of  the  work  at  UCLA  PPL  included 

•  set  up  and  operate  the  GPIM  dye  laser  in  a  wavelength  regime  suitable  for  probing 
the  (2-0)  band  of  the  first  positive  system  (750-800  nm), 

•  determine  the  precise  wavelengths  of  the  rotational  lines  of  this  band, 

•  determine  the  radiative  lifetimes  of  the  laser  excited  states, 

•  determine  the  linewidths  of  the  fluorescence  transtions, 

•  determine  the  photoabsorption  cross  sections  (using  the  above  measurements), 

•  re-evaluate  first  positive  LIF  lidar  performance  with  these  results,  and 

•  evaluate  the  performance  of  a  first  positive  LIF  lidar  utilizing  the  same  dye  laser  in 
the  GPIM  mobile  lidar  system  at  the  HIPAS  observatory. 

B.2  Description  of  Experiment 

Figure  4  shows  the  experimental  setup.  The  GPIM  dye  laser  is  pumped  with  a  UCLA 
PPL  XeCL  excimer  laser.  The  laser  output  is  routed  to  the  plasma  chamber,  located  in  a 
different  room  in  the  facility.  The  total  optical  path  is  approximately  35  m.  The  light  is 
directed  through  a  window  at  the  center  of  the  top  of  the  cylindrically  shaped  chamber. 
The  light  travels  through  the  chamber  along  the  diameter  of  the  cylinder,  through  another 
window  at  the  bottom,  and  hits  a  photodiode.  A  spectrometer  is  located  on  the  side  of 
the  chamber  so  as  to  detect  fluorescence  radiation  at  90°  scattering  angle  from  the  optic 
axis. 

The  dye  laser  is  a  Lambda  Physik  model  FL2002E.  This  laser  contains  an  oscillator  and 
preamplifier  sections  which  use  one  dye  cuvette,  and  a  main  amplifier  section  which  uses 
a  second  cuvette.  We  operated  this  laser  with  a  variety  of  dyes  during  the  course  of  the 


24 


Figure  4:  Experimental  setup  of  laser  induced  fluorescence  experiments  on  UCLA  PPL 
Racetrack  rf  plasma  chamber. 
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experiment  in  order  to  find  the  one  that  provided  for  the  most  energy  in  the  wavelength 
regime  of  interest.  The  dye  which  we  found  most  suitable  was  LDS765  (Exciton).  With  this 
dye  in  both  the  oscillator /preamplifier  and  amplifier  cuvettes  we  obtained  2.5-3  mJ/pulse 
in  the  range  760-780  nm.  We  obtained  1-2.5  mJ/pulse  in  the  range  740-760  nm.  These 
energies  were  obtained  without  the  use  of  the  manufacturer-supplied  intracavity  etalon 
in  the  oscillator.  With  the  etalon  the  energies  were  approximately  10%  of  these  values. 
According  to  the  manufacturer’s  specifications,  the  etalon  reduces  the  laser  linewidth  from 
0.2  cm-1  to  0.04  cm-1.  As  a  comparison,  the  linewidth  of  a  Nj  first  positive  rotational  line 
at  200° K  is  0.025  cm-1. 

The  excimer  laser  used  to  pump  the  dye  laser  is  a  Lambda  Physik  model  EMG-201 
operating  with  XeCl  at  308  nm.  The  output  energy,  which  produced  the  dye  laser  energies 
cited  above,  was  75  mJ/pulse.  The  pulse  repetition  rate  was  6  Hz. 

Most  of  the  experiments  were  performed  on  nitrogen  plasmas  produced  in  the  UCLA 
PPL  Racetrack  Experiment.  Nitrogen  plasmas  are  created  by  rf  breakdown  of  Nj  gas.  The 
gas  is  backfilled  into  the  chamber  to  pressures  in  the  range  of  10-4  to  10-s  Torr.  Different 
pressures  produce  different  plasma  conditions.  Most  of  the  experiments  were  performed 
with  10~s  Torr  backfill.  Approximately  10  kW  of  rf  energy  at  500  kHz  is  delivered  by  an 
antenna  located  inside  the  chamber.  The  rf  is  pulsed  at  6  Hz,  in  synchronization  with  the 
excimer  laser.  Some  experiments  were  performed  in  the  UCLA  PPL  Torch  experiment.  In 
this  experiment  Nj  can  be  backfilled  into  the  chamber  to  as  much  as  1  Torr,  producing 
much  higher  density  plasmas  than  in  Racetrack. 

The  spectrometer  used  to  detect  fluorescence  radiation  is  a  Jarrell  Ash  model  82-000. 
This  is  a  1/2  meter,  f/8  Ebert  type  spectrometer  (uses  a  grating).  The  minimum  resolvable 
wavelength  difference  in  the  740-780  nm  range  is  approximately  0.5  A. 

The  detector  was  a  photomultiplier  tube  operating  in  current  mode.  The  output  of  the 
PMT  was  directed  to  a  Tektronix  digitial  oscilloscope  with  200  MHz  maximum  digitization 
rate  and  averaging  capability. 

B.3  Experimental  Technique 

The  first  major  task  after  setting  up  the  dye  laser  to  operate  suitably  in  the  (2-0)  band  of 
the  first  positive  system,  was  to  observe  laser  induced  fluorescence.  Since  the  wavelengths 
were  calculated  theoretically  a  detailed  wavelength  scan  of  the  laser  was  required. 

Our  technique  for  observing  laser  induced  fluorescence  followed  the  following  steps: 
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1.  Tune  the  dye  laser  wavelength  through  the  (2-0)  band,  from  740  nm  to  790  nm, 
in  steps  of  0.003  A.  We  expect  this  step  size  to  be  more  than  adequate.  The  laser 
linewidth  is  approximately  0.2  cm-1,  or  0.1  A  at  7737  A,  without  the  etalon.  The 
linewidth  of  the  rotational  lines  of  the  neutral  molecules  in  the  experiment  was  about 
0.03  cm-1,  or  0.02  A  at  7737  A  (doppler  spread  at  room  temperature).  The  separation 
between  rotational  lines  in  the  first  positive  system  is  on  the  order  of  several  tenths 
of  A. 

2.  Observe  (wall)  scattered  laser  radiation  with  the  spectrometer  set  for  widest  pass 
band,  ~  10  A.  Doing  this  assures  us  that  the  spectrometer  is  set  for  the  same  wave¬ 
length  as  the  laser.  The  signal  from  the  scattered  radiation  is  approximately  10  nsec 
in  duration,  the  laser  pulse  width. 

3.  Shift  the  spectrometer  center  wavelength  as  the  laser  wavelength  is  scanned. 

4.  Observe  the  fluorescence  signal  due  to  the  decay  of  the  molecule  via  the  reverse 
pathway  of  the  laser  excited  transition.  The  signal  from  fluorescence  should  add 
an  exponentially  decaying  tail  to  the  signal  from  wall  scattered  radiation.  The  tail 
should  have  a  time  constant  equal  to  the  radiative  lifetime  of  the  upper  state  of  the 
fluorescence  transition  (~  10  nsec). 

5.  If  the  previous  step  is  successful,  then  with  the  same  excitation  wavelength  observe 
fluorescence  in  other  branches  of  the  same  band  and  in  other  bands  (e.g.,  the  (2-1) 
band  which  has  band  head  at  A  ~  8695  A). 

6.  Continue  with  the  wavelength  scan  to  map  out  the  rotational  structure  of  the  (2-0) 
band. 

With  the  above  data  lifetimes,  linewidths,  temperatures,  and  absorption  cross  sections 
can  be  determined.  Lifetimes  and  linewidths  come  directly  from  the  fluorescence  signals. 
Temperatures  are  determined  by  either  calculating  from  the  linewidths  of  individual  ro¬ 
tational  lines  or  by  considering  the  relative  strengths  of  signals  from  different  rototional 
lines.  The  latter  gives  the  rotational  temperature.  Cross  sections  are  determined  using 
Eq.  3,  with  transition  probability  (inverse  of  lifetime)  and  linewidth  from  the  data. 
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B.4  Results  and  Conclusions 


As  mentioned  above,  we  were  not  successful  in  demonstrating  laser  induced  fluorescence 
in  the  first  positive  system.  We  do  not  have  a  definitive  explanation  for  this,  although  we 
cite  some  possible  reasons. 

•  There  may  be  less  JVj(A)  than  we  expect  in  the  Racetrack  experiment.  This,  of 
course,  should  have  no  bearing  on  experiments  performed  in  the  ionosphere.  We 
estimate  N3(ASE+)  density  in  the  aurora  from  direct  observations  (see  Section  4). 

•  The  background  fluorescence  may  have  been  too  high.  If  this  is  the  case  then  it  is  * 
necessary  to  average  over  more  laser  pulses.  We  averaged  over  as  many  as  twenty 
pulses. 

•  The  absorption  cross  section  is  lower  than  we  expect.  This  is  the  only  reason  we  cite 
which  would  also  affect  our  lidar  calculations  in  this  report. 

We  have  some  evidence  that  suggests  the  first  reason  cited  above  may  be  the  real 
reason.  We  compare  expected  spectrometer  signals  from  second  positive  LIF  and  first 
positive  LIF.  Second  positive  LIF  was  observed  in  the  Racetrack  experiment.  The  ratio  of 
spectrometer  signals  is  given  by 

Pnn  l iVx.o\  ( Ra,j\  ( l (iy/ir)mr' \  / / rsc  \  ( Fb-°a\  ( Eqa\ 

Pssn  \Nb, 0/  \Rb,j)  V-Cfic(A)/  V(>?iC)s37l/  W/  VT.AB/ 

(Na,q\  (  01  \  (0.1  \  (0.1\  ( l-43  x  10~U\  (90  X  10"®\  /0.1617N  /3\ 

\Nb,o )  V0.05/  V  1  /  V0.3/  V  3.7  x  lO"18  )  \  9  x  10-«  J  V0.4527/  U ) 

=  X4.1  X  10"6 

\nb,oJ 

where  we  have  substituted  probable  values  for  each  of  the  quantities. 

If  the  first  positive  signal  to  be  equal  to  the  second  positive  signal  then  Racetrack  must 
have  2.4x10s  times  as  many  Ni(AsE+)  as  N](BsIIf).  This  may  not  be  the  case  in  the 
experiment.  In  any  case,  if  there  are  only  a  small  number  of  N3(j48E+)  that  simply  implies 
a  small  signal  in  which  case  more  sensitivity  is  required  of  the  detector/data  acquisition 
system  (e.g.,  averaging  over  more  shots). 

We  know  that  the  laser  was  scanning  the  correct  wavelength  regime.  Figure  5  shows 
passive  emission  from  the  plasma  as  a  function  of  wavelength.  The  data  was  acquired  with 
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the  spectrometer  near  the  end  of  the  rf  pulse.  The  spectrometer  was  set  for  approximately 
1  A  resolution.  Data  was  averaged  over  four  shots.  Different  bands  are  clearly  discernable. 
However  individual  rotational  lines  are  not  resolved. 

Appendix  C  Auroral  Emission  Data  Acquired  at  HIPAS, 

Alaska 

During  the  autumn  of  1992  and  winter  of  1993  passive  emission  from  aurora  was  observed 
using  the  receiver  of  the  Phillips  Laboratory  (Geophysics  Directorate,  GPIM  branch)  mo¬ 
bile  lidar  system  at  HIPAS  observatory  in  Alaska.  The  overall  purpose  of  these  observations 
was  to  evaluate  the  lidar  receiver  performance  in  possible  LI?  lid?r  wavelength  regimes 
(e.g.,  7737  A,  4278  A,  3914  A).  The  specific  goal  of  these  observations  were 

•  to  acquire  data  with  the  lidar  receiver  during  aurora  at  different  wavelengths, 

•  correlate  with  measurements  from  calibrated  photometers  operating  at  3914  A, 

•  estimate  lidar  sensitivity  at  potential  LIF  wavelengths,  and 

•  re-evaluate  previous  LIF  lidar  calculations. 

To  date  several  data  sets  have  been  acquired.5  Figures  6  and  7  show  two  sets  of  data. 

The  first  figure  shows  emission  at  774  nm  and  the  second  figure  shows  emission  at  428  nm. 

The  bandpass  for  both  cases  is  10  nm.  Both  figure  show  plots  of,  as  a  function  of  time, 
the  number  of  photons  collected  per  second  (left  hand  side)  or  the  column  emission  rate  in 
Rayleighs  (right  hand  side).  The  data  is  cut  off  as  the  photon  count  rate  increases  because 
of  sunrise.  Both  plots  indicate  emission  activity  above  background,  although  much  lower 
than  even  a  weak  aurora  would  produce. 


5Data  was  acquired  by  Richard  Dickman  of  the  HIPAS  observatory  and  by  Ralph  Wuerker  of  UCLA 
PPL. 
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(A*)  !•■*!»  1IU 


Figure  5:  Passive  emission  as  a  faction  of  wavelength  in  the  first  positive  system  from  the 
UCLA  PPLA  Racetrack  plasma. 
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PL/GRM  lidar  receiver  data 
HFAS  Observatory,  Alaska 
25  January  1993 


150  g 


CWL-774  nm 
BW—9.7  nm 


minutes  after  05:58:50  UT 


Figure  6:  Passive  emission  at  774  nm  recorded  by  the  GPIM  lidar  receiver  at  HIPAS 
observatory. 
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photons  collected/sec 


PL/GPIM  lidar  receiver  data 
HIPAS  Observatory,  Alaska 
15  March  1993 


minutes  after  08:14:02  UT 

CWL-428  nm 
BW=X)  nm 


Figure  7:  Passive  emission  at  428  nm  recorded  by  the  GPIM  lidar  receiver  at  HIPAS 
observatory. 
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